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Abstract 



Cosmological solutions for a pressureless fluid in the Brans-Dicke theory exhibit asymptotical 
accelerated phase for some range of values of the parameter ui, interpolating a matter dominated phase 
and an inflationary phase. The effective gravitational coupling is negative. We test this model against 
the supernovae type la data. The fitting of the observational data is slightly better than that obtained 
from the ACDM model. We make some speculations on how to reconcile the negative gravitational 
coupling in large scale with the local tests. Some considerations on the structure formation problem 
in this context are also made. 

PACS number(s): 04.20.Cv., 04.20.Me 

1 Introduction 

The accumulation of supernovae type la data seems to confirm that the universe is today in an accelerated 
phase 012] ■ The first indications of this accelerated phase came from the analysis of about 50 supernovae. 
Today, the data are approaching 300 supernovae and those first conclusions, for the moment, are 

confirmed. The theoretical explanation of this surprising result is one of the most important challenge 
in cosmology. In order to obtain an accelerated phase, a repulsive gravitational effect must be obtained. 
Generally, this is achieved by introducing an exotic fluid which exhibits a negative pressure today, called 
dark energy (for a recent review, see reference [3]). The most natural candidate is a cosmological constant, 
since it is an inevitable consequence of considering quantum fields in a curved space-time. But, the 
cosmological constant proposal faces a major difficulty: the theoretical predicted value surpass the value 
induced by observation by 120 order of magnitudes. This discrepancy can be alleviated, but not solved, by 
considering, for example, supersymmetric theories. Moreover, the introducing of a cosmological constant 
today does not explain the coincidence problem, the fact that the universe begun to accelerated very 
recently, after the completion of formation of local structures like galaxies and clusters of galaxies. 

Many other alternatives to the cosmological constant have been presented in the literature. In the 
quintessence program, a self-interacting scalar field is considered |HJ[7]- The potential for the scalar field 
is such that initially the pressure is positive, in order to allow the formation of local structures, becoming 
latter negative, driving the acceleration of the universe. The coincidence problem may be addressed in this 
case. However, fine tuning is required. However, some more natural potentials have been determined using 
supergravity theories Other dark energy models appear in the literature: k-essence UJ, Chaplygin 
gas ^HE]> viscous fluid J^l, etc. Here, we explore another possibility: the acceleration phase could be 
generated by a non-minimal coupling between a scalar field and the gravitational term. 

The paradigm of a scalar-tensor theory with non-minimal coupling is the Brans-Dicke theory, which 
is characterized by the Lagrangian 
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where a; is a coupling parameter and L m is a matter Lagrangian. Local tests at level of the solar system 
and stellar binary systems, restrict w to be larger than around 1,000 ^Sj- This makes the predictions 
of the Brans-Dicke theory almost indistinguishable of those coming from general relativity. The scalar 
field (f> is related to the inverse of the gravitational coupling G. Some recent interesting proposals address 
the possibility that the gravitational coupling can depend on the scale, having different value at local 
and at cosmological scales ^J^l- If such scale-dependent gravitational coupling can be satisfactorily 
implemented, it is possible to have smaller values for the parameter uj at cosmological scales, obtaining 
scenarios that differ substantially from those obtained from general relativity while agreement with local 
tests is preserved. Studies of quantum effects in gravitational theories, and their consequence for the 
behaviour of the gravitational coupling, have strengthened those proposals ^3 E| • 

In reference , general cosmological solutions for a flat universe in the Brans-Dicke theory were de- 
termined. For a pressureless fluid, such solutions can interpolate a decelerating phase with an accelerated 
phase. This is an example of how cosmological models in Brans-Dicke theory may differ substantially 
from the standard cosmological model: using general relativity such interpolation can just be obtained 
through introduction of an exotic fluid. However, the price to pay is not negligible: the Brans-Dicke 
parameter must takes values in the interval —3/2 < uj < —4/3, implying that the effective gravitational 
coupling is negative. The possibility of a scale-dependent gravitational coupling may render, however, 
such scenarios attractive. 

These solutions with a late accelerated phase were studied in reference ^Hl- Here we make a step 
further: a comparison with the supernovae type la data is made, using the "gold sample". This Brans- 
Dicke late accelerated model gives better agreement with the observational data than the ACDM model. 
The best fitting for the Brans-Dicke parameter and the Hubble parameter h are obtained, indicating 
uj w —1.5 and h w 0.6. The approach here differs from others existing in the literature which studies 
accelerated models in Brans-Dicke theory [23121112211231211231211123 by the fact that no other ingredient 
is introduced, like a potential term for the scalar field. But, the model requires an effective mechanism 
that leads to a repulsive gravitational coupling at cosmological scales, while keeping it attractive at local 
scales. 

This article is organized as follows. In next section, the pressureless Brans-Dicke cosmological model 
is revised. In section 3, the comparison with the supernovae type la "gold sample" is made. In section 4 
we present our conclusions. 



2 Late accelerated phase in the Brans-Dicke cosmological model 

Let us consider a flat Brans-Dicke cosmology, with the universe filled with pressureless matter. The 
equations of motion resulting from the Lagrangian Q are: 

»(:)' ■ t + S(?H? ■ < 2 > 

- i& • < 3 > 

P + 3-P = . (4) 
a 

The last equation leads to p = poa~ 3 . Inserting this relation in ©, we find the first integral, = J^fjC, 
where C is a constant. 

Following |28j . we can define an auxiliary function u, satisfying the relation 

± = -^+ 2 --" (5) 

u a t t w 

which, in view of the equations of motion, results in the integral relation 

du 



t[u + 4 ± s\/u 2 + Au] 



= 2In(t-t c ) , (6) 
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where s = 3y 1 + |w. This integral relation has three critical points u = 0, u = —4 and u = 4+ 2 3 ^ . The 

first one gives non physical results, while the second one leads to a oc t 2 for any value of w. For this last 
case, = —4(3 + 2uj). However, the fact that the evolution of the scale factor does not depend on lo 
makes this solution not very interesting. 

The third critical point of © leads to a particular solution in terms of power law function (choosing 
t c = 0): 

a = aut 1 ^ , 4> = (/> £3+3^ . (7) 

A more general solution may be obtained through integration of ©. This has been done in [18] . resulting 
in the general flat solutions 



a = ao(*-t-) 4 + 3 " (t-t+) , (8) 

</> = 0o(*-*-) 4+3 " (*-*+) 4+3 ^ . (9) 

where i± = t c + A(s ± 1), A being another integration constant. Since i+ > t_ , t = t+ may be identified 
with the initial time. 

The solutions presented above have the following properties. For t — > t+, they reduce to 



M , (io) 



l± 3v /l + l^ 

= <M , (11) 

while for t >> t + , they take the form 

a = a a ti^ , 4> = ^o* 31 ^" ■ (12) 
Concerning the possible scenarios, the solutions have the following properties: 

1. For uj > -|: The universe has a decelerated expansion during all its evolution. 

2. For lo < — | < uj < — | there are two regimes: 

(a) For the positive upper sign, initially the universe has a subluminal expansion, followed by a 
superluminal expansion; 

(b) For the negative lower sign, the universe has an accelerated phase during all its evolution. 

Is the case (2. a) that it will interest us. In such a case, the interpolation between a decelerated phase 
and an accelerated phase can be achieved. However, there is high price to pay in order to obtain such a 
simple qualitative realization of the accelerated expansion today without introducing any kind of exotic 
matter: the gravitational coupling is negative. In fact, the gravitational coupling is related to the scalar 
field (j> by the relation [2Hj 

G=—\ (13) 
4 + 3w (j) { ' 

Hence, for -3/2 < u < -4/3, G < 0. 

The most popular model to explain the current acceleration of the universe is the AC DM model. The 
matter content is a pressureless fluid and the cosmological constant. For a spatial flat case, the Einstein's 
equation reduce to 

/ • \ 2 

where p m = p m o& 3 - This equation admits the solution 



3 - =8ttGK + A , (14) 



/3A/V , 2 3A 
Kia) Smh *V4M* • (15) 
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where M = 8irGp m o/3. The asymptotical behaviour are 

t^O => aocti , (16) 

t — > oo a oc e v 4M . (17) 

As expected, initially the scale factor describes a dust dominated universe, and later a cosmological 
constant dominated universe. 

3 Comparing with the supernovae type la data 

We intend now to compare the background model for the flat pressureless fluid Brans-Dicke cosmology 
with the supernovae type la data. This is made through the computation of the luminosity distance 
function, given by 

D L = %, (18) 
a 

where do is the value of the scale factor today, a is the value of the scale factor at the time of the 
emission of the radiation, and r is the co- moving radial distance. From now on, we fix ao — 1. With this 
normalization choice, the scale factor at a given moment t is related with the redshift z by 

' (19) 



1 + 2! 

Hence, the luminosity distance can be expressed as 

D L = (l + z)r . (20) 
Usually, the observational data are expressed in terms of the distance moduli given by 

/i0 = 51 ° g (^) +25 ' (21) 
In general the task now would consist in using the Friedmann equation as 

> , (22) 



2 8nG^ 



where pi are the different matter components. If each of this component obeys a separate conservation 
equation, it can express as pi = pioa~ 3 ( 1+ai \ where a.i is the barotropic equation of state (supposed to 
be constant) defined by pi = oupi. Using then the relation between the scale factor and z, the luminosity 
function can be expressed as an integral over z which depends on the density parameter today for each 
fluid. This procedure is detailed, for example, in reference It must be remarked also that the 

expression (|21|) must be modified for the case the gravitational coupling varies with z, as reported in 
references [H01I21E3- However, in the case we consider here, the local physics is supposed to be given 
by general relativity, the Brans-Dicke modifications intervening at cosmological scales only. Hence, with 
this assumption, the expression l|21(l may be still used safely. 

For the Brans-Dicke cosmological model presented in the previous section the procedure outlined 
above becomes impossible to be applied due to the non-minimal coupling between gravity and the scalar 
field. But, we can use the relation 

a(t) 

where t e is the emission time for the source at the redshift z and to is the present time, together with the 
normalization condition a(to) = ao = 1. The expression JHJ may be rewritten as 

fx — — \ r+s 

o(x)= — tH x r ~ s , (24) 

^ to ' 



1 + §W t 

— , x = — . (25) 

3w 4 + 3w t+ 
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The present time may be determined through the definition for the Hubble parameter 

H= a^ 2rt - (r - s)t+ 

a (t-t+)t ' 

Evaluated for today, t = to, when H = Hq, this leads to the expression 



(H t+ + 2r) ± + 2r)2 - 4ff (r - s)t+ 
l °- 2H ' (2?) 

where Hq = 100 h km/ (M pes), h being one of the parameters to be determined. Hence, the luminosity 
distance is given by, 

D L = (1 + z)ct [ , (28) 

Jx a{y) 

where y is the integration variable. The same procedure is used for the AC DM model described in the 
end of the previous section. 

In order to compare with the supernovae data, the (dimcnsionless) time x of emission must be deter- 
mined from the expression for the scale factor and the relation between the scale factor and the redshift. 
In the case of the AC DM model, there are two parameters: h and A (or equivalently, p m o)- In the 
Brans-Dicke model, there are in principle three free parameters: cj, t + and h. However, we first inspect 
the best value for f + , and then we vary the parameters u> and h around this parameter. The best fitting 
are obtained for t+ = 5.2 x 10 18 : other values leads to an universe that accelerates too early or too 
late. Hence, in both models, we are left with two free parameters. We use the supernovae type la "gold 
sample", with 157 high redshift supernovae, described in reference 0] (see also reference |29p. 




Fi gure 1: Two-dimensional probabil- Figure 2: Probability distribution for Figure 3: Probability distribution for 
ity distribution for u and h. the parameter u. the parameter h. 




Figure 5: Probability distribution for Fig ure 6: Probability distribution for 
n m in the ACDM model. h in the kCDM model. 

Figure 4: Probability distribution for 
f2 m and h in the KCDM model. 

The main quantity to be evaluated is the \ 2 function, which gives the fitting quality: 



X 2 =J2 — < ( 29 ) 

i=i Gi 
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where /x*o is the predicted theoretical value for the ith supernova, /x° its observational value and u.- L the 
error bar, taking already into account the effect of peculiar dispersion. From this quantity, a probability 
function is obtained: 

P(u,h) = Ae-£ , (30) 

where A is a normalization constant. The probability function for a unique parameter can be obtained 
by marginalizing on the remaining one: 

p fP(u,h)dh JP(m,h)du J 

" J J P(iU,h)dhdLU ' H J J P(iU,h)dhduJ ' [ 1 

In figures 1, 2 and 3 the two dimensional probability distribution for lu and h, as well as the one- 
dimensional probability distribution for u> and h, in the Brans-Dicke flat model, are displayed. The 
most probable value for lu is —1.5, exactly the limiting case for the validity of solutions (JSJ. However, 
lu = —1.5 represents the case where the Brans-Dicke theory is conformally equivalent to the Einstein's 
theory. Hence, the prediction should more properly stated by saying that the most probable value for 
the Brans-Dicke parameter is u> — > —1.5. For h, the most probable value is h = 0.605. The best fitting 
for the Brans-Dicke flat model is given by lo — -1.477 and h = 0.594, with x 2 = 1.19318. 

In figures 4, 5 and 6 the two-dimensional probability for £l m and h and the one-dimensional probability 
function Q m and h, in the case of the ACDM model, are displayed. For h, the prediction is quite similar 
to the Brans-Dicke case: the most probable value, marginalizing over f2 m , is ft. = 0.597. The most 
probable value for the dark matter density parameter, marginalizing over h, is Q m — 0.436. The best 
fitting for the ACDM model is given by Sl m = 0.506 and h = 0.597 with X 2 = 1-19813. The Brans-Dicke 
flat model gives a fitting slightly (but not negligible) better than the ACDM model. 



4 Conclusions 

The general solutions for the Brans-Dicke flat cosmological model, determined by Gurevich et al, predict 
a late time accelerated universe for —3/2 < lu < —4/3, with an initial decelerated phase. This solution 
can be a candidate to describe the observed universe. In the present work, we have constrained this 
Brans-Dicke flat model using the supernovae type la " gold sample" . We have compared the results with 
those obtained in a simplified ACDM model. The Brans-Dicke model leads to a slightly smaller \ 2 (which 
characterizes the quality of the fitting of the observational data) with respect to the ACDM model. 

In both models, the predicted value for the Hubble parameter h is around 0.6. This contrasts strongly 
with the estimation of h coming from the anisotropy of the cosmic microwave background radiation, 
which leads to h = 0.72 ± 0.05 However, this seems to be a general feature of the estimation of h 
using supernovae data and CMB data [2SI • For the Brans-Dicke model, the best value for the parameter 
lu is around —1.5. It must be remarked that the case lu = —1.5 implies that the Brans-Dicke theory is 
conformally equivalent to general relativity. The minimum value for x 2 is obtained for h = 0.594 and 
lu = -1.477, with x 2 = 1-19318. In the ACDM case, the best fitting is obtained with f2 m = 0.506 and 
h = 0.597, with x 2 = 1-19813 

The results reported here indicate that, from the point of view of the supernovae type la data, the 
Brans-Dicke cosmological flat model leads to a quite viable scenario, even if the generalization to the non- 
flat cases should be made in order to have a more complete statistical analysis. However, at this point, 
this model must be seen as essentially a toy model, mainly due to the negative value for the parameter 
lu. It would be interesting to connect the Brans-Dicke model exploited here with effective actions in four 
dimensions coming from M-theory and F-theory, which predict also a negative value for the parameter 
to. 

The main problem with this scenario is that the values of lu are largely outside the values obtained with 
local tests, which indicate lu > 1,000 ^Sj- In reference (20j, indications for a negative lu have also been 
obtained. Moreover, the range —3/2 < lu < —4/3 implies a negative gravitational coupling j2H]. These 
drawbacks can be, in principle, circumvented if the gravitational coupling is scale-dependent, as suggested 
by considering quantum effects (161117) . In the present work, we made no attempt to reconcile the values 
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for uj deduced from local tests with the corresponding values deduced from the analysis performed here. 
But, evidently, this problem must be addressed in order to have a complete realistic scenario. 

One important point to be signed is that previous study indicate that structures can form in the 
Brans-Dicke model considered here during all the evolution of the universe, after the radiative phase, 
even the gravitational coupling is, at large scale, repulsive | 34l I35| . However, a more detailed comparison 
between the theoretical predictions for matter agglomeration and the observational data must be made 
in order to constraint more strongly the model. In special, the possibility that the model can be unstable 
due to the repulsive gravitational coupling must be verified 36 . 
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